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Abstract 

T-H ■ Cosmological reconstruction of Little Rip model in f{R,T) gravity is investigated, where R is 

■ the curvature scalar and T the trace of the energy momentum tensor. The model perfectly repro- 
^ : duces the present stage of the universe, characterized by the AC DM model, without singularity at 
^-H ' future finite-time (without the Big Rip). The input parameters are determined according to Su- 

' pernovae Cosmology data and perfectly fit with the WMAP around the Little Rip. Moreover, the 

o: 

_ thermodynamics is considered in this Little Rip cosmology and it is illustrated that the second law 

" ^ ^ I of thermodynamics is always satisfied around the Little Rip universe for the temperature inside the 

^ I 

5^ horizon being the same as that of the apparent horizon. 

X ■ 

U ; Pacs numbers: 04.50.Kd, 95.36.-^x, 98.80.-k 

1 Introduction 

Recent astronomical data from Type la supernovae 1 as well as from the CMB spectrum f2l confirm 
that our universe is undergoing an accelerated expansion period. In order to comply with this feature, 
dark energy content for the universe may be postulated |3], with undesired properties, such as the 
violation of some energy conditions. Other way is to consider extended theories of gravity, which besides 
their ability to describe the inflation in the early universe [4], can be used as candidates to explain the 
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present acceleration of the universe [5] without the requirement to introduce the exotic matter with 
negative pressure, dubbed as dark energy. A method to generahze the Einstein theory is to consider 
the gravitational Lagrangian as a function of the curvature R and the trace T of the energy momentum 
tensor, the so-called f{R,T) theory [6l. This theory of modified gravity has attracted more attentions 
and some interesting results have been found [7]-|llj. 

Recently, a novel scenario has been proposed in [291 , the so-called "Little Rip" (LR) . Models without a 
future singularity in which the energy density of the dark energy, pde, increases with time will nonetheless 
eventually lead to a dissolution of bound structures at some point in the future. In the LR, the scale 
factor and the density are never infinite at a finite time. Typically, the LR interpolates between the Big 
Rip (BR), where the scale factor and energy density diverge, and the ACDM model, where there is no 
such divergence and no disintegration because the dark energy remains constant. Mathematically, the 
LR can be represented as an infinite limite sequence which has the BR and the ACDM as its boundaries. 

In this paper, we consider the process of occurrence of the LR and propose to investigate the cor- 
responding explicit form of type R + 2g{T) model. Here, we are considering f{R,T) as f{R,T) = 
fi{R) + h{T), where we specially consider fi{R) = R and /2(T) = 2g(r). 

Note that black hole thermodynamics indicates the fundamental connection between the gravitation 
and thermodynamics |121 113| . The Einstein equation was derived from the Clausius relation in ther- 
modynamics with the proportionality of the entropy to the horizon area in General Relativity (GR) 
[l5] . This procedure was developed to more general extended gravitational theories [16l [17] , and also 
in so-called /(T) gravity [T3], where T is the torsion scalar. Note also that a new method is now 
developed for studying the thermodynamics of equilibrium system, the so-called geometrodynamics, and 
some interesting results have been found [T8]-|21j. 

In this paper the proposal is the same but in the framework oi R + 2g{T) gravity, analysing the validity 
of the law of thermodynamics in a LR model. In this optic, we investigate the non-equilibrium description 
of thermodynamics near the LR. We show that the second law of thermodynamics is always verified in 
LR model for the temperature inside the horizon being the same as that of the apparent horizon. 

The paper is organized as follows. In Sec. 2, we present the general formalism of f{R,T) gravity; 
the special R + 2g{T) gravity is adopted and the model leading to LR is obtained. In Sec. 3, the 
thermodynamics in LR cosmology is considered, where the first and second law are investigated. The 
conclusion is presented in the Sec. 4. 
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2 General formalism of f(R, T) gravity 

Let us consider the general model of f{R,T) gravity whose action can be described by the action 



(1) 



where k = 8ttG, with G the Newtonian gravitational constant, R is the curvature scalar and T the trace 
of the energy momentum tensor, defined from the matter Lagrangian density Cm by 

2 5{^C , 

Varying the action S with respect to the metric g'"^, one obtains [B], 



T,. - '\ir' ■ (2) 



fniR, T)R^, - -f{R, T)g^, + [g^.D - V^V,) fniR, T) - kT^, - MR, T)T^, ~ fM T)e^, , (3) 



where Qf^i, is defined by 

e - o"/3 ^'^"^ - 2T +a C 2q"/3 ^^^"^ (A) 

Here fji and fx denote the derivatives of / with respect to R and T, respectively. Let us assume for 
simplicity that the function / is given by f{R,T) — fi{R) + f2{T), where fi{R) and f2{T) are arbitrary 
functions of R and T, respectively. Then, ([3]) can be written as 

fm{R)R^.. - lfi{R)9^. + (ffM-n - V^V,) fMR) = ^T^. ~ f2T{T)T^. - hT{T)Q^, + i/2(r)g^, . (5) 

Assuming that the matter content is a perfect fluid, the stress tensor is given by 

T^v = {p + p)Uf,u„ - pg^,y , (6) 

where is the four- velocity satisfying u^u^ — 1. Then, the matter Lagrangian density can be taken as 
Cm = —p, and Q^i, — —2T^^ — pg^iv Hence, Eq. ^ becomes 



fiR{R)R^.. - ^/i(i?).9p. + (g^^n - v^v,) fMR) = kT^, + /2T(r)T^, 



/2T(r)p + i/2(T) 



ffM^ ■ (7) 



Now, we propose to assume a special case where fi{R) = R and f2{T) = 2g{T), i.e., f{R, T) = R + 2g{T) 
as 2g{T) additive term to Einstcin-Hilbcrt one. This case seems interesting and has been widely studied 
for other purposes [5]-|10) and interesting results have been obtained. Then, Eq. ([7]) becomes 

i?^!. - ^Rg,,u = kT),^ + [2gT{T)p + g{T)] , k = k + 2gT{T) . (8) 

We see from this equation that coupling constant k in GR, becomes a running constant k in i? + 2g{T) 
gravity. 

Let us now consider the spatially flat Friedmann-Robertson- Walker (FRW) line element 

ds^ ^dt^ -a'^{t)dK^ , (9) 
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where a{t) is the scale factor. One easily gets the first and second generalized Friedmann equations, 
respectively as 



-2H - iH^ = k 



1 



p + - C^pgr + g) 

K 

P - ^ (2j3.gT + g) 



(10) 
(11) 



where H = a/a is the Hubble parameter and the dote denoting the derivative with respect to the cosmic 
time t. Calling the expressions in brackets in the right side of ^TU\\ and ([TT|) respectively the effective 
energy density pef / and pressure p^f / , one write ((T0| and (fTTj) as 



3H'^ = kpeff , 



-2H - 



KPeff 



(12) 
(13) 



Note that by combining ([T^ and ([T^. one can write peff = —Peff — 2H/k. Since, the effective energy 
density is directly linked with the square of the Hubble through ((T2|) . one can take H as function of Peff- 
Then, one can write the equation of state as 



Peff 



-Peff - fiPeff) 



(14) 



From the conservation law, the effective energy density evolves as 

d{kpeff) 



dt 



= -3Hk{p^ff +peff) 

= 3Hkf{peff) . 



Let us now assume that the scale factor can be written as [27] 



a{t) 



(15) 



(16) 



where h{t) is a non-singular function. Hence, the effective energy density can be written as Pe// — ih^ /k. 
The condition for kp^ff being an increasing function of the scale factor a is d{kpeff)/da — {6/d)hh > 0. 
This holds as long as h > 0. By using Eq. (|15p. it appears that an increasing kp^ff is ensured by 
kf{peff) > 0. Integrating Eq. ((Tsj) . one gets 

r""" d{kp) 
Ipo 3k/(p)^ 



a = ao cxp 



(17) 



Also, Eq. ([H]) yields 



t = 



d{kp) 



Piff 



(18) 



The condition for a big rip singularity is that the integral converges. Choosing / such that f{p) 
A{p/kf/^ EH [28], one gets 



Inf^ 



(19) 
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A simple analysis of ()19|) shows that npeff does not diverge until an infinite time t has elapsed. This is 
the Little Rip phenomenon. From (|17|) . we can express Pe// in functions of a, 



Peff = -jPo 



1 



and also using Eq. (|12p . one obtains after integration 



a(t) = oo exp ■ 



iA 



exp — - — t — 1 



(20) 



(21) 



that we put in a compact form as 

a{t) = aoexpck (e''* - l) , (22) 
where a and (3 can be easily found by identifying ([22|) with (|2T|) . as a = 2^Kpo/(3yl) and (3 = A^/3/2. 

2.1 Little Rip + 251 (T) type model 

Here, we propose to construct model of type R + 2g{T) according to Little Rip scale factor ((22|) . 
With ([m, the Hubble parameter and its first derivative read 



H{t) = a/3e^* , H = afi^e^' . 
Using the relation p = poa~^^^~^^^ and P^ . the trace is written as 

T = p - 3p = po(l - 3a;) exp [-3a(l + a;)(e'^* - 1)] 



(23) 



(24) 



where the equation of state p = up is assumed for the ordinary content of the universe. Substituting the 
second term of ((23|) into ((24|) . one acquires 

H = + \n\[pfi{\-Zuj)]^^) TWTz:)~^ . (25) 

On the other hand, by summing (|T0| with pT|) . one gets 

p(l+a;) (k + 2gT) + 2ij = 

By using and the first equality of ([M)) . Eq. becomes 
1 + w 



(26) 



1 - 3a; 

whose general solution reads 



{k + 2gT) + 'iaP^ + 2\n I [po(l - 3a;)] TaiT+iio 



(27) 



.(T) = -- 



T a/32(l-3a;) 



k(1 +0;) 



ln<^ [po(l - 3a;)] ^5^1+1:7 T3(i+i:o 



3(1-80;) 
2^ 



In^ <^ [po(l - 3a;)]^TT+:^ I + C- 



(28) 



where Ci is an integration constant. The corresponding R + '2g{T) model is 
m T) = R-T- ^"^y - In ( [pod - 3.)] ^ 

K{l+UJj [ 

+ ^^^^^^|^ln'|[po(l-3c^)]^TW^| +2Ci . (29) 

The constant Ci can be determined in the following way. For GR without cosmological constant, one 
has SHq = KpQ, where Hq and po are respectively the current values of the Hubble parameter and energy 
density. From the observational data, the current value of the Hubble parameter is Hq = 2.1 x 0.7 x 
10^'^^GeV [311I3S]. In this general case, working with f{R,T) model, the initial condition has to be the 
same, that is 3Hq = kpef f\t=to- This leads to 

2po(l+w)gT(To)+g(ro) =0 , (30) 

where Tq = Pq{1 — Slu) is the current trace of the energy momentum tensor. By using ([50]) and (|28p . one 

gets 

"1 - 3cl; + k(1 +0;) 



To . (31) 



2(1 - Scj) 

Furthermore, by using the initial condition npa = 3Hq, one obtains a = 2Ho/{AV3) and /3 = A%/3/2, 
which leads to a/3 — Hq. In [37], the parameter A, which also appears here is Eq. (|19l) is chosen to 
make a best fit to the lasted supernova data [22], and can be found for the range bounded as 2.74 x 10^'^ 
Gyr^^ < A < 9.67 x lO^^Gyr^^. We also have the constant n = SttG = 8tt/Mpi, with the Plank mass 
of AIpi ^ G"i/2 = 1.2 X lO^^GeV. Then, the parameters a and ^ can perfectly be found according to 
observational data. 

By calculating the parameter of effective equation of state, one obtains 

~Kp + 2pgTiT)+g{T) 



^eff 



kp + 2pgT{T)+g{T) 



kp + 2pgT{T)+g{T) 

= (32) 
3a 

It is easy to observe from l\'32\i that for any t from this present time to to the future, tOeff < — 1. This 
means that the universe is always in phantom phase. One can also observe that in the limit of large t, 
{t — > oo), uJeff — 1 and hence the Little Rip scenario can hold. 

An interesting feature to be noted in this model is that, from for small value of the cosmic 

time, meaning that we are around the present stage of the universqj, one gets e*^* — 1 ^ f3t and then, 



^The small values of time are values closed to the present time to, since, t may be substituted by t — to such that for t 
closed to to, the difference t — to is small. 
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a{t) cx e°('=''*-i) e"*^* = e"°K In this situation, H = and ujeff = -1; this is the ACDM model, where 
the gravitational action reads i? — 2A, with A the cosmological constant. 

In general this R+2g{T) may be view as a type of the so-called A{t) gravity where A{t) = 2g[T{t)], 
and this is an interesting model because it provides the interaction between matter and dark energy |24| . 
This model has an advantage over f{R) gravity theories because we do not need to transform (conformally) 
from a metric to another, and there is no question about which metric is physical since we use one and 
the same metric tensor |25| . 

Without loss of generality, let us assume that the Hubble parameter can be written as 

H{t) = HLRe-^^'---'^ , (33) 
from which the scale factor can be written as 



a{t) = aLR exp 



Hlr f^.p(t^^-t) _ ^ 



(34) 



where H^r and a^R denote respectively the Hubble parameter and the scale factor at Little Rip time 
tbR- Note that in the LR model, even if H is finite in finite future, in the same, when H becomes infinite, 
the inertial force becomes very strong and destroys any bounded object [28) . This has been dubbed as 
LR [57]. For the model of type p3l) . the time t^R — to from the present until the destruction of the 
Earth-Sun system has been estimated to t^R ~ to ^ 146Gyrs [28| . 

Then, one has the model that describes very well the present accelerating universe, and presents any 
finite time singularity in future, i.e., the Big Rip. 



3 Thermodynamics in LR cosmology 

Here, we present thermodynamics in f{R, T) and try to check the viability of such a theory. By 
following the same procedure in Refs. |30[ 131] . we examine whether the second law of thermodynamics 
can be verified near the LR. In GR, the proportionality of the entropy to the horizon area allows to 
derive the Einstein equation from the Clausius relation in thermodynamics |32| . Consequence of this has 
been applied to more general extended gravitational theories |33[ I34j , and we propose to perform this in 
f{R,T) gravity around the Little Rip. 

3.1 First Law of Thermodynamics 

According to the recent type la Supernovae data, it is suggested that in the accelerating universe the 
enveloping surface should be the apparent horizon rather than the event one from the thermodynamic 
point of view ^35| . In flat FRW spacetime, the radius at the apparent horizon is given by = 1/H [36] , 
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and the dynamical apparent horizon is determined by the relation 

h^^dafdpf = , Kfi = diag (l, -a^it)) . (35) 
The first derivative of the f a , leads to 

dfA = -r\HHdt . (36) 
On the other hand, the equations ^T2\\ and ([T3|) can be put into the following set of equations 

- f Pe// , (37) 

H = -^(Peff+Peff) ■ (38) 

By making use of ([38|) . one can write (p6|) as 

dfA = ^f\H {p,ff + p,ff) dt . (39) 

In GR, the Berkenstein-Hawking horizon entropy is expressed as Sbh — ■Aj (4G) [Hl[35], where A = 47rr^ 
is the area of the apparent horizon. Note from Eq. ([8]) that the Newtonian gravitational constant G is 
replaced by the running gravitational coupling parameter G = G + 2gT{T)/{8T:). Then, in the context of 
this special modified gravity, as in [32], the entropy is defined as Wald entropy and reads S — A/{4:G). 
By using pop. the derivative of the horizon entropy is written as 

dS - 87r\\H ip^ff + p^ff) dt-^dG . (40) 

One can now calculate the Hawking temperature = |(Tsg |/(27r), corresponding to the associated 
temperature of the apparent horizon, where asg denotes the surface gravity which is given by |39] 

1 / rA _^ 



TA 



= -^[l-iu:,fj)p,fj , (41) 

where we used Vef f I Pef f — ^eff- We clearly see from (|4T1) that the surface gravity is positive for 
uJeff > 1/3. Hence, the temperature is obtained as 

Th = ^(i-^) . (42) 



2TrrA V 2i/fA 

Now, by combining (|42|) and (j40|) . we get 



TndS = iirr^H {p.fj + p^//) dt - 2^f\H {p.jf + p.fj) d^A + "^d^A ~ ^dG . (43) 



Defining th.e Misner-Stiarp energy as 



S = ^ - yPeff , (44) 



where V — Anf\/i. 

From the equation of continuity, one acquires d(pe//)/di + 3i/ (pe// + Pe//) = —3H'^k/{K^). By using 
this, one gets 

dE = A7rr\Hp,ffdfA - ^7Tf\H (p.ff + p,ff) dt - ^dG . (45) 
Defining the work density as [40], W = -{l/2)T°'f^'^^f ^'^h ap, one acquires 

W=l{Peff-Peff) ■ (46) 

Now, using dV = Airf'^drA, one can perform the product WdV as 

WdV = 2nP\{p,ff-p,ff)dPA ■ (47) 



By summing (|43| with p5)) . one easily obtains (|47|) plus an additive term, that is TndS = —dE + WdV — 
TndSad, which can be written as 

TudS + TadSad = -dE + WdV , (48) 
TndSad = ~^ ^ '^^^^ dG . 

This is the non-equilibrium description of the thermodynamics. This additive term may be interpreted 
as an entropy production term in the non-equilibrium R + 2g{T) gravity. Near the LR, this additive term 
should be view as the entropy linked with the dissolution of bound structures in LR. This result is quite 
similar to that obtained in [41]|43]. Observe that if g{T) — 0, the gravitational action is reduced to the 
Einstein-Hilbert's one and the equilibrium relation is achieved. Specially in this paper, we propose to 
analyse the behaviour of the entropy in LR model. One can remark from (|40l) that 

^2-3 



S = Sn r^{peff +Peff) 



Q2 



2PaH + £ ) . (49) 



"'A 

G \ " G 



By using Eq. (gH]), Eq. (gSl, Eq. and the relation 8ttG ^ k = k + 2gT, one gets 

G^-^-0-^^.Me'^^Ml^^){e^^~l)^ , (50) 

which means that the running gravitational parameter G is negative in this LR model. Let us show 
clearly the main reason for which G is negative. Observe that by using Eq. one gets 



In GR the gravitational action term is just the Einstein-Hilbert term R, if the cosmological is note 
considered. In such a situation, in order to obtain an accelerating expanded universe, one needs to 
consider an exotic component (the dark energy), with negative pressure such that the parameter of 
equation of state is negative {lode < 0). However, when the gravitational action is modified, one acquires 
the dark energy effect from the contribution of the modified gravity terms, in our case, 2g{T). Since 
this contribution has to play the role of dark energy, just the parameter ujdf associated to the co-called 
dark fiuid may be negative (cudf < — !)• One has to consider that the weak energy condition is always 
satisfied for the ordinary matter, that is p + p >, or 1 + oj > 0. Moreover, the LR model is guaranteed 
for H ^ h > 0. By using these evidences, it follows that the quantity k + 2gT is negative, which means 
that G < 0. 

Moreover, by using ([50]). one acquires 



G 
G 



[/3 + 3a(l + w)] e'^* > 



(52) 



Note that H — h, and the condition for the occurrence of the LR is /i > 0. Thus, from (|49| . one sees 
that in the LR universe and the framework of i? + 2g{T) gravity, the horizon entropy S always increases 
S>0. 

In GR where g{T) — 0, one has G — G ^ const and the expression (|49| reduces to 

-2nH 



Sgr — 



GH^ 



(53) 



In this case, for phantom universe (H < 0), the entropy always decreases. This also corresponds to the 
situation where the null energy condition is violated, i.e. pde + Pde < 0. This is also occurs in some 
type of modified gravity |36| . In modified gravity the sign of the first derivative of the horizon entropy 
depends on the form of /, i.e. the form of f{R,T), or f{R) or f{G), depending of what type of modified 
gravity is used {G is the Gauss-Bonnet invariant term) [46l |47] . 

One can now express the horizon entropy corresponding the R + 2g{T) model compatible with the 
LR, as 



S 



A 

4G 



am 

Sn'^Poil+uj) 



exp [-3l3t - 3a(l + w) (e'^* - l)] 



Sn^il +uj)T 



(3a; - 1) [a/3 + k{T)] [aP"^ + kiT)] 



k{T) = In 



(54) 
(55) 



It is very important to note here that in the framework of GR, the LR cosmology is such that the horizon 
entropy and its first derivative read 

TT TT ofl, • 27r 



S 



-2/3t 



cm Ga^P^ 



GPa- 



:e-2^* < 



(56) 
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Once again, we see from this that the behaviour of the entropy quite depend on the form of the grav- 
itational action. Hence, we see that, stiU guaranteeing the accelerated expansion of the universe, the 
behaviour of the entropy may change in a LR cosmology, depending of the form of the gravitational 
action in consideration. However, in any case, as the time goes to infinity, the entropy vanishes. 

Without any loss of generality, the assumption ([33]) may be assumed. One can write the horizon 
entropy around the LR as 



Slr = Sq exp 



-3/3(tLi^-to)-3«(l+L^)(e^(*-«-*°)-l)] , (57) 



By using ([551) . one can rewrite the entropy at the LR as 

r / Hr r.\ 11 

(59) 



5 = 5'oexp<j -31n I ] -3a(l+w) 



Ho 



ln(^)-l 



Slr = So { ^) = Soe'f^^'+'^^'+-^^ . (60) 



By considering that the LR must occur very far from today, one can use Hlr >> Ho and then, (|59p 
becomes 

-3[l+a(l+w)] 

Ho , 

Then, since we know the values of the parameters a and /3, the horizon entropy at LR may be calculated. 
Particularly, if we consider that the ordinary content of the universe is dust, i.e. w = 0, the horizon 
entropy at LR reads 

SLR{dust) = - nr,TT ^ ■ (61) 

3.2 Second law of thermodynamics 

In this subsection we consider the Gibbs equation in terms of all matter and energy fluid to examine 
the second law of thermodynamics, that is 

TndSin = d{peffV) +PeffdV = Vdpeff + (peff + Peff) dV , (62) 

where Sin denotes the entropy of all the matter and energy inside the horizon. By using the equation of 
continuity, ()62|) can be written as 

THdS^n = -AtT^aH {peff + Peff) + 47rf^ (p^ff + Peff) df A " ^^JG . (63) 

By summing (l63t with (|43)) . one gets 

Th {dS + dSin) = ^TTTA [Peff + Peff) df A + ^f^df^ - ^^'^ ' (64) 
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from which, by dividing by dt, and using the relation dvA/dt = k (peff ~^ Peff) one gets obtains 

2GIP HG^ 



The second law is described by 



dS dSi n „ , N 



Since the temperature Th is always positive, one gets the following condition for the second law of 
thermodynamics 

-^<- (77 + 77) • (67) 



H \H H 



As we have previously shown, the quantity G/G is always positive in this LR model, and both H and H 
are also positive. Then, ((67|) always hold for the model R + 2g{T), meaning that the second is always 
satisfied near the LR. By choosing the LR time « ^lr ^^'^ the present time as ~ ^cT^' easy 
to observe that —0.21 < e^'3(*i«^*o)/Q! < 0.165, in perfect agreement with the WMAP data ^2,, that is 
— 1.11 < ujde < —0.86. This proves the consistency of the result. 



4 Conclusion 

In this paper we reconstructed the cosmological model of effective phantom type and which does not 
lead to a singularity at finite future time (Little Rip model) in the framework of a type of /(i?, T) gravity, 
where R is the curvature and T the trace of the energy momentum tensor. Physically, in the LR, the 
scale factor and the energy density are never infinite at finite time. The model is reconstructed by solving 
an inhomogeneous differential equation of first order from the combination of the generalized Friedmann 
equations in the special case, R + 2g{T), where g{T) is only function of T. The integration constant 
Ci of the general solution is determined by imposing the same initial condition in GR to this modified 
gravity. Since the initial value of the trace To can be directly related to the current Hubble parameter 
Hq = 2.1 X 0.7 X 10^^^ GeV, the LR models is perfectly reconstructed. By the use of the Supernova 
Cosmology Project observational data, we show that the parameters a and /3 can be found from the know 
range of values of the parameter A, and then, the input parameters are perfectly in agreement with the 
observational data. Furthermore it shown that the model reproduce the ACDM model at this present 
stage of the universe and confirm that the LR interpolates between the ACDM model and the Big Rip. 

Moreover, we undertake the thermodynamics in this LR model in the framework of the special R + 
2g{T) gravity showing the non-equilibrium description of the thermodynamics. We also observe that the 
second law of thermodynamics is always satisfied in this model when the temperature inside the horizon 
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is the same as that on the horizon. An interesting point to be mentioned is that at LR time the current 
value oi cjde in this LR model is consistent with the WMAP observational data. 

Acknowledgement: M. J. S. Houndjo thanks CNPq-FAPES for financial support. M. E. Rodrigues 
thanks UFES for the hospitality during the elaboration of this work. 



References 

[1] T. M. Davis et al, Astrophys. J. 666, 716 (2007) |arXiv:astro-ph/0701510] . 

[2] Dunkley et al [WMAP CoUaboration] , Observations: arXiv: 0803.0586 [astro-ph], E. Komatsu et 
all. [WMAP Collaboration], arXiv: 0803.0547 [astro-ph]. 



[3] T. Padmanabhan, AIP Conf. Proc. 861, 179 (2006) [arXiv: |astro-ph/0603114| . 



[4] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980); M. B. Mijic, M. S. Morris and W. Suen, Phys. Rev. 
D 34, 2934 (1986); P. A. Anderson and W. Suen, Phys. Rev. D 35, 2940 (1987); T. P. Sotiriou, and 
V. Faraoni, Rev. Mod. Phys. 82, 451 (2010). 

[5] A. G. Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter et al, Nature (London) 391, 51 (1998); 
P. M. Garnavich et al, Ap. J. 509, 74 (1998); S. Permutter et al, Astrophys. J. 517, 565 (1999). 

[6] T. Harko, F. S. N. Lobo, S. Nojiri and S. D. Odintsov, '7(i?,r) gravity," Phys. Rev. D84 (2011) 
024020. (arXiv:1104.2669l [gr-ac]]. 

[7] M. J. S. Houndjo, " Reconstruction of f(R, T) gravity describing matter dominated and accelerated 
phases". Int. J. Mod. Phys. D. 21, 1250003 (2012). arXiv: 1107.3887 [astro-ph.CO]. 

[8] M. J. S. Houndjo and O. F. Piattella, "Reconstructing f(R, T) gravity from holographic dark energy" , 
Int. J. Mod. Phys. D. 21, 1250024 (2012). arXiv: 1111.4275 [gr.qc]. 

[9] D. Momeni, M. Jamil and R. Myrzakulov, Euro. Phys. J. C 72, arXiv: 1107.5807[physics.gen-ph]. 

[10] M. J. S. Houndjo, C. E. M. Batista, J. P. Campos and O. F. Piattella, "Finite-timne singularities in 
f{R,T) gravity and the effect of conformal anomaly". |arXiv: 1203.6084 [gr-qc]]. 

[11] F. G. Alvarenga, M. J. S. Houndjo, A. V. Monwanou and Jean. B. Chabi-Orou, "Testing some f(R, 
T) gravity models from energy conditions" arXiv: 1205.4678 [gr-qc]. 

[12] J. M. Bardccn, B. Carter and S. W. Hawking, Commun. Math. Phys. 31, 161 (1973); J. D. Beken- 
stein, Phys. Rev. D 7, 2333 (1973); S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum- 
ibid. 46, 206 (1976)]; G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2738 (1977). 



13 



[13] T. Padmanabhan,I^S^09l03l3p[gr-qc]; AIP Conf. Proc. 1241, 93 (2010) |arXiv:0911.i403l [gr- 
qc]]; Kept. Prog. Phys. 73, 046901 (2010). 

[14] M. Hamani Daouda, Manuel E. Rodrigues and M. J. S. Houndjo, Euro. Phys. J. C 72, 1890 (2012). 
arXiv: 1109.0528 [physics. gen-ph]. 

[15] T. Jacobson, Phys. Rev. Lett. 75, 1260 (1995). 

[16] E. Ehzalde and P. J. Silva, Phys. Rev. D 78, 061501 (2008); K. Bamba, C. Q. Ceng, S. Nojiri and 
S. D. Odintsov, Europhys. Lett. 89, 50003 (2010); S. F. Wu, B. Wang, X. H. Ge and G. H. Yang, 
Phys. Rev. D 81, 044010 (2010); Y. Yokokura, larXiv:1106.3149l [hep-th]. 

[17] R. Brustein and M. Hadad, Phys. Rev. Lett. 103, 101301 (2009). 

[18] Manuel E. Rodrigues and Zui A. A. Oporto, Phys. Rev. D 85, 104022 (2012). arXiv: 1201.5337[gr-qc]. 

[19] Manuel E. Rodrigues and Glauber T. Marques, "Thermodynamics of a class of non-asymptotically 
flat black holes in Einstein-Maxwell-Dilaton Theory", arXiv: 1206:0763[gr-qc]. 

[20] Deborah F. Jardim, Manuel E. Rodrigues and Stephane J. M. Houndjo, "Thermodynamics of phan- 
tom Reissner-Nordstrom-AdS black hole", arXiv: 1202.2830[gr.qc]. 

[21] Manuel E. Rodrigues, Deborah F. Jardim and Stephane J. M. Houndjo, "Thermodynamics of black 
plane solution," arXiv: 1205. 3481[gr-qc]. 

[22] R. AmanuUah, et aI,(The supernovae Cosmology Project) Ap. J 716, 712-738 (2010). 

[23] N. J. Poplawski, rarXiv:gr-qc/0608031| >r2. 



[24] M. Bronstein, Phys. Z. Sowjetunion 3, 73 (1933); M. S. Berman, Phys. Rev. D 43, 1075 (1991). 



[25] G. Magnano, gr-qc/9511027 (1995). 



[26] Makarenko Audrey N., Obukhov Valery V. and Kirnos Ilya V. " From Big to Little Rip in modified 
f(R, G) gravity", arXiv: 1201.4742 [gr-qc]. 

[27] Paul H. Frampton, Kevin J. Ludwick and Robert J. Scherrer, " The Little Rip", Phys. Rev. Rev. D 
84, 063003 (2011). arXiv: 1106. 4996 [astro-ph.CO]. 

[28] I. Brevik, E. Elizalde, S. Nojiri and S. D. Odintsov, " Viscous Little Rip". Phys. Rev. D 84, 103508 
(2011). arXiv: 1107.4642 [hep-th]. 

[29] P. H. Frampton, K. J. Ludwick and R. J. Sherrer, Phys. Rev. D 84, 063003 (2011). 

[30] K. Bamba and C. Q. Geng, JCAP 1111, 008 (2011) ; arXiv:1109.1694l [gr-ac]]. 



14 



[31] K. Bamba and C. Q. Ceng, JCAP 1006, 014 (2010) |arXiv: 1005.5231 [gr-qc]]. 
[32] T. Jacobson, Phys. Rev. Lett. 75, 1260 (1995) |arXiv:gr-qc/9504004j . 

[33] C. Eling, R. Guedens and T. Jacobson, Phys. Rev. Lett. 96, 121301 (2006) |arXiv:gr-qc/0602001] ; 
E. Elizalde and P. J. Silva, Phys. Rev. D 78, 061501 (2008) |arXiv:0804.372ll [hep-th]] : K. Bamba, C. 
Q. Ceng, S. Nojiri and S. D. Odintsov, Europhys. Lett. 89, 50003 (2010) [arXiv:0909.4397l [hep-th]]: 
S. F.Wu, B.Wang, X. H. Ge and G. H. Yang, Phys. Rev. D 81, 044010 (2010) |arXiv:0909. 13671 
[gr-qc]]; Y. Yokokura. larXiv:1106.3149l [hep-th] . 

[34] R. Brustein and M. Hadad, Phys. Rev. Lett. 103, 101301 (2009) arXiv:0903.0823l [hep-th]]; R. 
Brustein and A. J. M. Medved. iarXiv:1201. 57541 [hep-th]. 



[35; 
[36 
[37 
[38; 



[39 

[4o; 



[41 

[42; 
[43; 
[44; 
[45; 

[46 
[47 



J. Zhou, B. Wang, Y. Gong and E. Abdaha, Phys. Lett. B 652, 86 (2007). 

K. Bamba, R. Myrzakulov, S. Nojhi and S. D. Odintsov, l arXiv:1202.4057| i^3 [physics. gen-ph]. 

Brustein, R., Gorbonos, D. and Hadad, M.: Phys. Rev. 079(2009)044025. 

J. M. Bardeen, B. Carter and S. W. Hawking, Commun. Math. Phys. 31, 161 (1973); J. D. Beken- 
stein, Phys. Rev. D 7, 2333 (1973); S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum- 
ibid. 46, 206 (1976)]; G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2738 (1977). 

R. G. Cai and S. P. Kim, JHEP 0502, 050 (2005) |arXiv:hep-th/0501055] . 

S. A. Hayward, Class. Quant. Grav. 15, 3147 (1998) |arXiv:gr-qc/9710089] ; S. A. Hayward, S. 
Mukohyama and M. C. Ashworth, Phys. Lett. A 256, 347 (1999) [arXiv:gr-qc/9810006|; R. G. Cai 



and L. M. Cao, Phys. Rev. D 75, 064008 (2007) arXiv:gr-qc/0611071| 



Kazuharu Bamba and Chao-Qiang Ceng, larXiv;1109.1694v2 [gr-qc] . 
Brustein, R., Gorbonos, D. and Hadad, M.: Phys. Rev. D 79, 044025 (2009). 
M. Sharif and M. Zubair, [arXiv:1204.0848k ^2 [gr-qc]. 

E.W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, Redwood City, Cahfornia, 1990). 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) larXiyTTOOl 745381 
[astro- ph. CO]]. 

H. Mohseni Sadjadi, "On thermodynamics second law in the modified Gauss Bonnet gravity" , 
larXiv:1009.1839^2 [gr-qc] . 

H. Mohseni Sadjadi, Phys. Rev. D 73, 063525 (2006). 



15 



